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ARTICLE INFO ABSTRACT 


Keywords: 
Camellia sinensis 


Pu-erh tea has become one of the popular and widely consumed teas because of its health benefits resulted from 
various bioactive compounds that are microbially produced during Pu-erh tea production. We applied 16S rRNA 
Fungi gene amplicon sequencing and analytical tools to investigate microbiomes and associated chemical profiles of 


ces five raw and five ripened Pu-erh teas from different manufacturers. Contrasting bacterial and fungal micro- 
_ ae oe biomes of raw and ripened Pu-erh tea samples associated with distinct chemical profiles were revealed. The 
‘u-erh tea 


ripened Pu-erh teas had a higher bacterial diversity but lower fungal diversity than the raw Pu-erh teas. 
Firmicutes and Proteobacteria were the dominant phyla in the ripened and raw Pu-erh tea, respectively. 
Cyanobacteria phylum was found very abundant in the raw but rarely present in the ripened Pu-erh tea. 
Ascomycota was the dominant fungal phylum in both ripened and raw Pu-erh tea. Tea polysaccharide and 
gallocatechin gallate were the characteristic chemicals in the ripened Pu-erh tea; while caffeic acid, caffeine, 
EGC, ECG, EGCG, C, and GA were more representative in the raw Pu-erh tea. These results will assist in further 
cultivating and characterizing the representative microorganisms in Pu-erh teas of different origins, and improve 


our understanding of microbial-mediated Pu-erh tea production process. 


1. Introduction 


Pu-erh tea is a popular dark tea produced in southern Yunnan of 
China (Liang Zhang, Zhang, Zhou, Ling, & Wan, 2013), and known to 
benefit consumers due to its hypolipidemic, anti-obesity, anti-micro- 
bial, anti-mutagenic, anti-oxidative, anti-cancerogenic, free radical 
scavenging, and toxicity suppressing activities (Lee & Foo, 2013; Y.; 
Zhang et al., 2016). These beneficial health effects are related to various 
bioactive compounds (e.g., polyphenols, caffeine, and purine alkaloids), 
which are produced by microorganisms during the fermentation pro- 
cess (Liang Zhang et al., 2013). Depending on the aging and fermen- 
tation processes, Pu-erh tea can be classified into two types: raw and 
ripened. The raw Pu-erh tea is processed as a typical green tea, and 
subject to the natural fermentation process during the storage. The ri- 
pened Pu-erh tea differs from the raw Pu-erh tea in an additional “pile 
fermentation” step involved with water replenishment to shorten the 
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aging process, resulting in the distinct chemical and microbial profiles 
between raw and ripened Pu-erh teas (Y. Zhang et al., 2016). 

Pu-erh tea is produced from the microbial fermentation of dried and 
rolled leaves of tea plant Camellia sinensis var. assamica (Y. Zhang et al., 
2016). Diverse microorganisms (e.g., bacteria, fungi) dominate dif- 
ferent steps of Pu-erh tea production and affect the production quality 
of Pu-erh tea. Unravelling the Pu-erh microbiomes is critical for de- 
veloping safe, consistent and high-quality Pu-erh teas. 

Our understanding of the Pu-erh microbiome has greatly advanced 
with the development of molecular tools. Earlier studies relied on tra- 
ditional methods (e.g., isolation, cultivation) to investigate the micro- 
organisms in the Pu-erh tea (Q. Wang, Gong, Chisti, & 
Sirisansaneeyakul, 2015). Fungi were considered as dominant and in- 
dispensable microorganisms in Pu-erh microbiomes, and many fungal 
genera, such as Aspergillus and Saccharomyces, have been isolated from 
Pu-erh tea samples (Z. J. Zhao, Tong, Zhou, Wang, & Liu, 2010). 
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Nevertheless, the microbial communities of Pu-erh tea are still not 
completely deciphered due to inadequate samplings, inappropriate 
methods of identification, and incapability of isolating uncultivable 
fungi and bacteria (Z. J. Zhao et al., 2010). By combining traditional 
cultivation-based methods with cultural-independent methods, the 
microbial diversity contributing to Pu-erh tea has been revealed gra- 
dually. For instance, a total of 320 bacterial and 112 fungal isolates 
have been identified with both culture-dependent and -independent 
methods, and bacterial diversities of Pu-erh teas were found to be po- 
sitively correlated with aging time, while fungal diversity of Pu-erh teas 
peaked after the sixty months of aging (Tian, Zhu, Wu, Wang, & Liu, 
2013). 

Compared with traditional cultural-independent methods, metage- 
nomic-based high-throughput sequencing (e.g., amplicon sequencing) 
has the advantages of increased sampling depth, high-throughput in a 
single run, and cost-effective. Amplicon sequencing (e.g., 16S rRNA 
genes) was applied to investigate microbial communities of different 
fermented foods (e.g., cheeses, teas) (Aldrete-Tapia, Escobar-Ramirez, 
Tamplin, & Hernandez-Iturriaga, 2014; Chen, Chen, & Lei, 2017; De 
Filippis, Parente, & Ercolini, 2017; Delcenserie et al., 2014; Zhong 
et al., 2016). For example, an amplicon sequencing study revealed a 
total of 629 bacterial and 390 fungal operational taxonomic units 
(OTUs) in the tea leaves, raw, and ripened Pu-erh teas, and showed 
aging only changed the microbiomes of raw Pu-erh teas but not the 
ripened Pu-erh teas (Y. Zhang et al., 2016). Another lab-scale study 
with metagenomic sequencing identified 2319 bacterial and only 39 
fungal OTUs in the Pu-erh tea, and Pseudomonas and Aspergillus were 
dominant bacterial and fungal genera, respectively (M. Zhao et al., 
2015). 

Catechins (e.g., especially (—)-epigallocatechin gallate (EGCG), 
(—)-epicatechin gallate (ECG)), theanine and sugars contents decreased 
greatly because of microbial digestion during the Pu-erh tea fermen- 
tation process. Production of methoxybenzene derivatives was probably 
performed by Aspergillus niger, Aspergillus clavatus, and Basidiomycetes 
spp. (Gong et al., 1993). Although chemical constituents of Pu-erh tea 
have been characterized by many studies, new tools and _ inter- 
disciplinary studies continue bringing novel insights into the chemical 
profiles of Pu-erh tea. Additional efforts are required to untangle the 
complex chemical constituents of Pu-erh teas and the association be- 
tween the microbiomes and chemical profiles of raw and ripened Pu-erh 
teas, thus helping us better understand the production process of Pu-erh 
teas. 

This study is aimed at identifying the representative microorgan- 
isms in the raw and ripened Pu-erh tea, establishing an associated re- 
lationship between chemical and microbial profiles of Pu-erh tea, and 
promoting our understanding of microbial-mediated Pu-erh tea manu- 
facture processes by profiling and comparing the microbial commu- 
nities and chemical compositions in Pu-erh tea samples collected from 
different manufacturers in Yunnan province. 


2. Materials and method 
2.1. Sample collection 


Five raw (R1-R5) and five ripened (F1-F5) Pu-erh tea cakes were 
purchased from eight manufacturers in the Menghai region of Yunnan 
province. Triplicate samples were harvested and analyzed from each tea 
cake. The details about the tea cakes are listed in Table 1. The tea 
samples were collected from the factories and then transported to the 
laboratory in sterilized polyethylene bags. All samples were stored at 
— 20 °C until further processing. 


2.2. DNA extraction, library construction, and sequencing 


DNA extraction from the Pu-erh tea samples was performed fol- 
lowing the established protocols (J. Wang et al., 2019). Briefly, after 
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DNA extraction, the quality of DNA samples was evaluated by elec- 
trophoresis on 1% agarose gels. The V4 region of the 16S SSU rRNA 
gene and fungal internal transcribed spacer (ITS) region were targeted 
by the primer sets 515F-806R and ITS1F-ITS2R, respectively. Each DNA 
sample was amplified in triplicate PCR reactions with a total volume of 
20 uL. The PCR products were pooled together and purified with the 
AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, 
USA). Purified PCR products were quantified by QuantiFluor™ (Pro- 
mega, Madison, WI, USA) and mixed. Then an Illumina paired-end li- 
brary was constructed using the pooled DNA by following the Illumina 
genomic DNA library preparation procedure. The subsequent 250 bp 
paired-end sequencing was performed on the Illumina MiSeq platform 
(San Diego, CA, USA) at Shanghai BIOZERON Co., Ltd. according to the 
standard in-house protocols. The raw sequence reads were deposited 
into the NCBI Sequence Read Archive (SRA) database under the sub- 
mission number of SUB5897881. 


2.3. Sequencing processing and analysis 


The raw sequencing reads were processed with the in-house pipeline 
by using Trimmomatic (Bolger, Lohse, & Usadel, 2014), FLASH (Magoc 
& Salzberg, 2011), Usearch (Edgar, 2010), and Qiime (v1.7.0) 
(Caporaso et al., 2010) software for denoising, trimming, quality fil- 
tering, pairing, and aligning. Usearch software (version 7.1, dri- 
ve5.com/uparse) was used to cluster the operational taxonomic units 
(OTUs) at 97% similarity, and then the representative sequence of each 
OTU was assigned to taxonomic information using the RDP classifier 
with a 70% threshold by comparing with Silva database (Release119, 
www.arb-silva.de) for prokaryotes and Unite database (Release 6.0, 
www.unite.ut.ee) for fungi. Community richness index (e.g., Chao1), 
community diversity (e.g., Shannon Index, Simpson Index), and the 
Good's coverage of sequencing were calculated using the Mothur Soft- 
ware (Schloss et al., 2009). Differences in the microbial community 
structures among Pu-erh tea samples were assessed using phylogeny- 
based weighted UniFrac distance or taxonomy-based Bray-Curtis dis- 
tance. Principal coordinate analysis (PCoA) was performed using the 
“ade4” package (Dray & Dufour, 2007) of R software (R_Core_Team, 
2017). Linear discriminant analysis (LDA) of effect size (LEfSe) was 
applied to determine the most discriminant taxa among Pu-erh tea 
samples of different manufacturers. The “anosim” module in “vegan” 
package of R software was used to conduct anosim analysis (Oksanen 
et al., 2017). Taxonomic and phylogenetic tree of microbial flora in Pu- 
erh tea was presented using the GraPhlAn software (Asnicar, Weingart, 
Tickle, Huttenhower, & Segata, 2015). 


2.4, Chemical compound analysis 


One gram of tea powder (40 mesh) was steeped in hot water (90 °C) 
for 40 min, and followed by filtration. The extraction was repeated 
three times, and the supernatant was mixed with 80% ethanol. After 
4 h’ precipitation, the mixture was filtered, and the residue was dis- 
solved with hot water and diluted to 100 mL for further use. Tea 
polysaccharides were quantitated using the anthrone-sulfuric acid 
method (Laurentin & Edwards, 2003). 

Tea polyphenol was determined by the ferrous tartrate method 
(Sava, Yang, Hong, Yang, & Huang, 2001). The contents of theanine 
(TA), caffeine (CAF), chlorogenic acid (CGA), gallic acid (GA), caffeine 
(CAF), caffeic acid (CA) and catechins, including (+)-catechin (C), 
(—)-epicatechin (EC), (—)-epigallocatechin (EGC), ECG, (—)-Galloca- 
techin gallate (GCG), (—)-Gallocatechin (GC) and EGCG in the Pu-erh 
tea samples were determined by high-performance liquid chromato- 
graphy (HPLC) using an Agilent 1200 series HPLC System with a diode 
array detector (DAD) (Santa Clara, CA, USA) following the published 
method (He et al., 2015). All chemicals standards (Purity =>98%) were 
purchased from Yuanye Biological Technology Co., Ltd. (Shanghai, 
China). Two gram of tea samples were extracted with 50 mL 50% 
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Table 1 
Pu-erh tea sample information. 
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Tea Cake Type of Pu-erh Tea Date of Manufacture Tea Factory Address 
R1-1; R1-2; R1-3 Raw 2016/09/19 Yunnan Yipintang Tea Industrial Co.,Ltd. (#1) Menghai County, 
Yunnan, China 
R2-1; R2-2; R2-3 Raw 2017/02/20 Menghai Tea Factory (#2) Menghai County, 
Yunnan, China 
R3-1; R3-2; R3-3 Raw 2017/05/18 Menghai Yinyun Tea Factory (#3) Menghai County, 
Yunnan, China 
R4-1; R4-2; R4-3 Raw 2017/03/08 Menghai Taixiang Tea Factory (#4) Menghai County, 
Yunnan, China 
R5-1; R5-2; R5-3 Raw 2017/03/18 Menghai Fuhai Tea Factory (#5) Menghai County, 
Yunnan, China 
F1-1; F1-2; F1-3 Ripened 2017/08/06 Putnam Pu-erh tea processing factory (#6) Ninger County, Pu-erh City, Yunnan, China 
F2-1; F2-2; F2-3 Ripened 2017/02/22 Menghai Tea Factory (#2) Menghai County, 
Yunnan, China 
F3-1; F3-2; F3-3 Ripened 2017/08/21 Menghai Huapu Tea Industrial Co.,Ltd. (#7) Menghai County, 
Yunnan, China 
F4-1; F4-2; F4-3 Ripened 2017/01/08 Lao Zhai Tea Processing Factory (#8) Wuhua District, Kunming City, 
Yunnan, China 
F5-1; F5-2; F5-3 Ripened 2017/03/18 Menghai Fuhai Tea Factory (#5) Menghai County, Yunnan, China 


ethanol solution at 80 °C for 20 min. After suction filtration and freeze- 
drying, the extracts were re-dissolved in 35% methanol solution for 
HPLC analysis. Before injection in HPLC, samples were diluted fourfold 
with methanol and filtrated through a 0.45 um Millipore filter. The 
separation was carried out using an Agilent Zorbax SB-C18 column 
(5 um, 4.6 mm x 250 mm, Santa Clara, CA, USA). The mobile phase 
consisted of phase A: methanol and phase B: methanol-0.1% TFA so- 
lution. Linear gradient elution was run at a flow rate of 0.8 mL/min. 


2.5. Statistical analysis 


Statistical analysis was performed by JMP Pro version 13.2 (SAS 
Institute Inc, Cary, North Carolina, USA). ANOVA was carried out to 
determine significant difference (P < 0.01). 


3. Results 
3.1. Sequences statistics 


The total numbers of qualified bacterial raw sequences for the raw 
and ripened Pu-erh tea samples were 692,054 and 666,119, respec- 
tively (Table S1). By comparison, 709,040 and 680,205 qualified fungal 
sequences were obtained for the raw and ripened Pu-erh tea samples. 
The average numbers of qualified sequences eligible for downstream 
analysis of bacterial and fungal community were 45,272 + 9181 and 
46,308 + 6408 (Table S1), which are more than the numbers of se- 
quences reported by previous studies (Ma et al., 2017; M.; Zhao et al., 
2015). 


3.2. Microbial richness and biodiversity 


The quality-filtered bacterial and fungal sequences of all thirty Pu- 
erh tea samples were clustered into 1461 bacterial OTUs and 836 fungal 
OTUs at the 97-similarity cut-off, respectively (Fig. S1). The raw and 
ripened Pu-erh tea samples shared 763 bacterial OTUs and 98 fungal 
OTUs (Fig. S1). Although the average number of bacterial OTUs in 
fifteen ripened Pu-erh tea samples was slightly larger than that in fif- 
teen raw samples (Fig. 1A), the average bacterial OTUs numbers be- 
tween raw and ripened samples had no statistically significant differ- 
ence (P = 0.44). By comparison, more fungal OUTs were found in raw 
Pu-erh tea samples (Fig. 1E, P = 0.017). The Chaol index suggests a 
higher bacterial richness (Fig. 1B, P = 0.004) in the ripened Pu-erh tea 
samples compared with the raw Pu-erh tea samples. There was no 
statistically significant difference in the fungal richness between raw 


and ripened Pu-erh tea samples (Fig. 1F, P = 0.116). Both Shannon and 
Simpson indices indicate that the ripened Pu-erh tea samples had a 
higher bacterial diversity but lower fungal diversity than the raw Pu- 
erh tea samples (Fig. 1C, D, G, H). The fungal rarefaction curves 
reached the saturation plateau (Fig. S2); while the bacterial rarefaction 
curves did not fully approach the saturation plateau (Fig. S3). The 
Good's coverage estimates for the bacterial and fungal diversity in all 
Pu-erh samples were higher than 99%, suggesting that a major part of 
the diversity in all samples had been captured. 

To compare the bacterial and fungal communities between raw and 
ripened Pu-erh tea samples, the phylogeny-based metric, weighted 
UniFrac distance was used to assess beta-diversity. The Principal 
Coordinate Analysis (PCoA) reveals that the bacterial (Fig. 2A) and 
fungal (Fig. 2B) communities were distinctly different between the raw 
and ripened Pu-erh tea samples, indicating that the manufacturing 
processes may have decisive impacts on the microbial community 
structures of the Pu-erh tea samples. Two outliers of raw Pu-erh tea 
samples (e.g., R5-3 and R1-1 samples in Fig. 2A) were relatively close to 
the ripened Pu-erh tea samples, probably due to the unexpectedly high 
abundances of Firmicute in R5-3 and R1-1 raw Pu-erh tea samples 
(Fig. 3A). 


3.3. Microbial community structures 


3.3.1. Bacterial community comparison 

The predominant bacteria phyla differed between the ripened and 
raw Pu-erh tea samples (Fig. 3A). Firmicutes was much more abundant 
in the ripened Pu-erh tea samples (58.2%) than the raw Pu-erh tea 
samples (13.1%) (P = 2.6 x 10~7); while Proteobacteria was the 
dominant phylum in the raw Pu-erh samples (60.53%) compared with 
ripened Pu-erh raw samples (14.7%) (P = 2.1 x 107°). Cyanobacteria 
phylum (19.37%) was found very abundant in the raw Pu-erh tea 
samples, while Cyanobacteria was rarely present in the ripened Pu-erh 
tea samples. At the genus level, Staphylococcus (22.7%), Lactococcus 
(15.9%), Bacillus (11.9%), Brevibacterium (7.3%), Saccharopolyspora 
(4.0%), Nesterenkonia (3.9%), Vibrio (3.8%), Kocuria (3.5%), Brachy- 
bacterium (2.6%), Pantoea (2.4%), Ralstonia (1.7%), Microbacterium 
(1.2%), Pseudomonas (1.2%), and Streptococcus (1.2%) were frequently 
detected in the ripened Pu-erh tea samples (Fig. 4A). In contrast, the 
dominant genera in the raw Pu-erh tea samples were Ralstonia (20.5%), 
Chloroplast norank (19.3%), Pseudomonas (14.2%), Lactococcus (7.0%), 
Halomonas (4.4%), Vibrio (3.3%), Bosea (2.1%), Escherichia/Shigella 
(1.8%), Ruegeria (1.8%), Bacteroides (1.2%), Prevotella 9 (1.1%), and 
Cellulophaga (1.1%) (Fig. 4A). 
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Fig. 1. The number of operational taxonomic units (OTUs) and alpha diversity indices for bacterial (A, B, C, and D) and fungal (E, F, G, and H) communities in the 
raw (red) and ripened (blue) Pu-erh tea samples. The bacterial communities in the ripened and raw Pu-erh tea samples are comparable diverse; while higher fungal 
diversity in the raw Pu-erh tea samples than in the ripened Pu-erh tea samples. (For interpretation of the references to color in this figure legend, the reader is referred 


to the Web version of this article.) 


3.3.2. Fungal community comparison 

More than 99.9% of fungal sequences in the ripened Pu-erh tea 
samples were assigned to the Ascomycota phylum (Fig. 3B). Similarly, 
Ascomycota (83.9%) was the dominant fungal phylum in the raw Pu- 
erh tea samples, followed by Basidiomycota (7.1%), Mucoromycota 
(1.8%), and Chytridiomycota (0.03%) (Fig. 3B). The fungal commu- 
nities in the ripened Pu-erh tea samples were mainly consisted of six 
genera: Blastobotrys (58.5%), Aspergillus (16.4%), Rasamsonia (10.0%), 
Thermomyces (7.5%), Candida (5.0%) and Strigula (1.73%); while the 
raw Pu-erh tea samples harbored more diverse fungal genera, including 
Aspergillus (29.0%), Cladosporium (9.8%), unclassified (7.2%), Eurotium 
(5.0%), Fusarium (3.8%), Rhodotorula (3.6%), Acremonium (3.6%), Pe- 
nicillium (3.3%), Blastobotrys (1.7%), Mortierella (1.7%), Schizothecium 
(1.7%), Kazachstania (1.6%), Debaryomyces (1.5%), Setophoma (1.3%), 
and Pestalotiopsis (1.1%) (Fig. 4B). 


3.3.3. Factors affecting microbial community 

We further analyzed the raw and ripened Pu-erh samples from the 
same factories: R2 and F2 were from Menghai Tea Factory (MTF), and 
R5 and F5 were from Menghai Fuhai Tea Factory (MFTF). Both PCoA 
and H-cluster tree analyses suggested that manufacturing factories did 
not play a role in separating the bacterial and fungal microbiomes, but 
the manufacturing processes (raw vs ripened) determined the clustering 
of bacterial and fungal microbiomes (Fig. S4). Anosim analysis con- 
firmed that there was no difference between samples from MTF and 
MFTF (R statistic = 0.094, P = 0.194 for bacterial microbiome; R 
statistic = 0.1, P = 0.148 for fungal microbiome). LEfSe analysis fur- 
ther revealed the most discriminant bacteria in the Pu-erh samples from 
the MTF and MFTF factories (Fig. S5), for instance, Blautia, Dorea, 
Megamonas, Loktanella, Chryseobacterium, Prevotella, Coprococcus and 
Pantoea for MTF tea samples, and Pseudomonas, Lactococcus, 
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Fig. 2. Principal coordinate analysis (PCoA) of bacterial (A) and fungal (B) communities in the raw (blue squares) and ripened (red dots) Pu-erh tea samples. The raw 
and ripened Pu-erh tea samples possessed distinct bacterial and fungal communities, which separated the Pu-erh tea into two clusters according to the PCoA analysis. 
(For interpretation of the references to color in this figure legend, the reader is referred to the Web version of this article.) 


Aerosphaera, Mycobacterium, Peptostreptococcus, Caldicoprobacter, Ste- 
notrophomonas, Carnobacterium, Geobacillus, Helcococcus, Moheibacter 
and Acinetobacter for MFTF samples; while Tausonia was the only re- 
presenting fungal group for the MFTF samples. 


3.3.4. Beneficial and harmful microorganism detection frequency in Pu-erh 
tea samples 
Fifteen beneficial microorganisms, such as Lactobacillus, Aspergillus, 
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Lactococcus, Penicillium, were compared for their abundances in the raw 
and ripened Pu-erh tea samples. Beneficial microorganisms were found 
more abundant in the ripened Pu-erh tea samples compared with the 
raw Pu-erh tea samples (Fig. 6SA). Conditionally pathogenic micro- 
organisms, such as Candida, were generally detected in low abundances 
in the Pu-erh tea samples (Fig. 6SB). Surprisingly, Candida parapsilosis 
was found abundant in the raw samples R4-2, R4-3, R5-1, and R5-2. 
Comparatively, Candida blankie was frequently detected in the ripened 
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Fig. 3. Bacterial (A) and fungal (B) community of the raw and ripened Pu-erh tea samples at the phylum level. Firmicutes and Proteobacteria were the predominant 
bacterial phyla in the ripened and raw Pu-erh tea samples, respectively. Ascomycota was the major fungal phylum in both the ripened and raw Pu-erh tea samples. 
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Fig. 4. Heatmap of bacterial (A) and fungal (B) communities in the raw and ripened Pu-erh tea samples at the genus level. The bacterial and fungal communities were 


different at the genus level for the raw and ripened Pu-erh tea samples. 


Pu-erh tea samples. 
3.4. Chemical profiles of Pu-erh tea 


The chemical compounds, including tea polysaccharide (TPS), CAF, 
TA, CA, CGA, catechins (e.g., GC, EGCG, ECG, EGC, C, EC, GCG, GA), 
and TPS were analyzed for all raw and ripened Pu-erh tea samples 
(Table 2). The average concentrations of TPS, CAF, TA and polyphenols 
extracted from the raw Pu-erh tea samples were 7.68 + 0.31 mg/g, 
30.74 + 0.23 mg/g, 78.61 + 4.11 mg/g, and 269.23 + 6.06 mg/g. 
The mean concentrations of TPS, CAF, TA and polyphenols extracted 
from the ripened Pu-erh tea samples were 35.50 + 0.33 mg/g, 
20.21 + 0.55 mg/g, 18.29 + 0.41 mg/g, and 159.43 + 2.11 mg/g. 
The differences of CGA, EC, C, and EGCG contents between raw and 
ripened Pu-erh samples were 1.50 mg/g (P = 0.0109), 6.44 mg/g 
(P < 0.0001), 2.31 mg/g (P = 0.0026), and 11.37 mg/g (P = 0.0292), 
respectively. No significant differences were found between raw and 
ripened Pu-erh tea samples for CA, GA, GCG, EGC, ECG, and GC con- 
tents. The principal component analysis suggested that the raw and 
ripened Pu-erh tea samples were clustered into two separate groups 
(Fig. 5). TPS and GCG were the characteristic chemicals found in the 
ripened Pu-erh tea samples; while the other chemicals (e.g., CA, EGC, 
CAF, ECG, EGCG, C, CA, GA) were more representative in the raw Pu- 
erh tea samples. 


3.5. Correlation between chemical compounds and microbiomes 


We explored the correlations between chemical compounds (e.g., 
TPS, TA, CGA, CAF, and catechins) and microbiomes at the phylum 
level (Fig. 6) by using the predictor screening method. Actinobacteria 
played a relatively important role in predicting the concentrations of 
TA, CAF, and TPS, while Cyanobacteria for CGA and GA, Bacteroidetes 


for CA, Gemmatimonadetes for GCG, Ascomycota for EC and C, some 
unclassified bacteria for EGC, Spirochaetae for ECG, and Firmicutes for 
EGCG and GC, respectively (Table 3). 


4. Discussion 


Microbiomes of Pu-erh teas greatly affect chemical compositions of 
Pu-erh tea. Unravelling the microbial communities responsible for 
producing beneficially bioactive compounds is critical in Pu-erh tea 
manufacture. We conducted experiments to investigate the chemical 
and microbial profiles of raw and ripened Pu-erh teas, and found con- 
trasting microbiomes of raw and ripened Pu-erh teas associated with 
distinct chemical profiles. 

Comparison of Chemical Profile Analysis. In this study, con- 
centrations of TPS, CAF, TA, catechins, and polyphenols in the raw and 
ripened Pu-erh tea samples were analyzed. Contents of tea polyphenols 
and catechins (e.g., EGC, C, EC, EGCG, ECG) were lower in the ripened 
Pu-erh teas than in the raw leaves, and the concentration of GA in- 
creased significantly in the ripened Pu-erh tea samples, which are 
consistent with a previous study (M. Zhao et al., 2015). . Ku et al. found 
the distinct chemical compositions between raw and ripened Pu-erh tea 
samples using a metabolomics approach, and reported higher anti- 
oxidant activities due to increased contents of total phenolic and fla- 
vonoid compounds in raw Pu-erh teas than in ripened Pu-erh teas (Ku, 
Kim, Park, Liu, & Lee, 2010). Zhang et al. applied the HPLC-DAD-ESI- 
MS" system to analyze fourteen ripened and fourteen raw Pu-erh tea 
samples quantitatively, and discovered the levels of total catechins in 
raw Pu-erh tea samples (i.e., 38.69 + 16.51 mg/g) were much higher 
than those in ripened Pu-erh tea samples (i.e., 4.28 + 1.01 mg/g). 
Gallic acid content, by comparison, was much higher in the ripened Pu- 
erh tea samples (L. Zhang, Li, Ma, & Tu, 2011). The average con- 
centration of caffeine in seven Pu-erh tea samples was 37.1 mg/g, 
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which is comparable to our findings; but the mean catechins content 
quantified in Liang's study was only 8.41 mg/g, which is much lower 
than what we measured (Liang, Zhang, & Lu, 2005). Differences in the 
chemical concentrations and profiles of Pu-erh tea samples observed 
among previously discussed studies are probably resulted from distinct 
age, manufactures, production processes of Pu-erh tea, and non-stan- 
dardized analytical methods. For instance, some studies sampled the 
raw Pu-erh tea several years after the production, and comparatively, 
we sampled the raw Pu-erh teas immediate they were produced in the 
factories. Natural fermentation occurred during the prolonged storage, 
which caused the differences of the chemical and microbial composi- 
tions among different Pu-erh tea samples. 

High through-put Sequencing Revealing Details of Pu-erh Tea 
Microbiomes. Different properties between raw and ripened Pu-erh 
teas (e.g., textures, aromas, color, flavors, taste, nutrients, chemical 
profiles) were caused by the indigenous microbiomes and distinct tea 
manufacturing processes. The indispensable roles of diverse micro- 
organisms (e.g., bacteria, fungi) in the raw and ripened Pu-erh tea 
production processes have been gradually unveiled (Z. Li et al., 2018; 
Lyu et al., 2013; Ma et al., 2017; Y. Zhang et al., 2016; M. Zhao et al., 
2015). A culture-based study found the predominance of yeasts and 
Aspergillus in the raw and ripened Pu-erh tea and identified 62 fungal 
isolates, including 13 species of Aspergillus, and seven species of Peni- 
cillium (Z. J. Zhao et al., 2010). Aspergillus niger and Blastobotrys ade- 
ninivorans were indicated to affect the nutritional enhancement of tea 
leaves during Pu-erh tea fermentation (Abe et al., 2008; W.; Zhang 
et al., 2016). These earlier studies were limited by inadequate sampling, 
non-optimized microbial isolation and identification techniques, and 
lack of microbial characterization (Z. J. Zhao et al., 2010). High- 
throughput sequencing technique offers a great tool to understand the 
microbiomes responsible for the production of Pu-erh tea. Lyu applied 
the pyrosequencing technique to profiling the bacterial and fungal 
community changes during the production of ripened Pu-erh teas, and 
found significantly lower percentages of Eukaryota (about 16%) and 
dominance of Bacteria (about 75%) in the ripened Pu-erh teas; they 
further demonstrated the dominant bacterial groups were Proteo- 
bacteria, Actinobacteria, and Firmicutes, while the major fungal groups 
included Saccharomycetes (8.06%), Yarrowia (1.49%), Saccharomyces 
(1.00%) and Aspergillus (0.86%) (Lyu et al., 2013). Another sequencing 
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Fig. 5. Principal component analysis of chemical 
Se profiles used to distinguish the raw and ripened Pu- 
erh teas. The raw and ripened Pu-erh tea samples 
* were separated into two clusters by TPS, TA, CAF, 
\ and catechins. TPS and GCG were the characteristic 
\ chemicals found in the ripened Pu-erh tea samples; 
\ while the other chemicals were more representative 
in the raw Pu-erh tea samples. 


study found genus Aspergillus was the dominant fungi, accounting for 
94.98% of the total fungal sequences (M. Zhao et al., 2015). The 
number of bacterial OTUs ranging from 948 to 2307 and the number of 
fungal OTUs ranging from 30 to 162 were identified in twelve fer- 
mented Pu-erh tea samples, respectively (Ma et al., 2017). We also 
applied the high-throughput sequencing techniques to investigate and 
compare the microbial diversities of raw and ripened Pu-erh tea, and 
found that Aspergillus and Blastobotrys were the two most abundant 
fungal genera in the ripened Pu-erh tea samples, while the raw Pu-erh 
tea harbored more diverse fungal genera. 

Microbiomes Associated with Chemical Profiles of Pu-erh Teas. 
Aspergillus plays a critical role in the production of Pu-erh tea; for in- 
stance, Aspergillus was suggested to be the principal theabrownin-pro- 
ducing microorganism during the early fermentation process (Q. Wang, 
Gong, Chisti, & Sirisansaneeyakul, 2014; Q. Wang et al., 2015). Despite 
that Blastobotrys was found very abundant in all fifteen ripened Pu-erh 
tea samples of our study and also previous study (Z. Li et al., 2018), this 
genus was not linked to the flavor-producing function. It has been hy- 
pothesized that thermotolerance and halotolerance Blastobotrys plays 
an indispensable role in the biological transformation of the tea leaves 
by applying monooxygenases and dioxygenases to oxidizing poly- 
phenols in the tea leaves (Abe et al., 2008). The reason for the in- 
creasing abundance of Blastobotrys and decreasing abundance of As- 
pergillus in the ripened Pu-erh tea is that water content in the 
fermentation pile is continuously reduced during the fermentation, and 
the late-stage environments favor the growth of Blastobotrys but not 
Aspergillus (H. Li et al., 2018). Penicillium was found relatively abundant 
in the raw Pu-erh tea samples but not in the ripened Pu-erh tea samples. 
There was a positive correlation between the abundance of Penicillium 
to the contents of (—)-epigallocatechin 3-O-gallate, (—)-epicatechin 3- 
O-gallate and 1,4,6-tri-O-galloyl-B-p-glucose (Ma et al., 2017). It is 
possible that Penicillium could secrete some antibacterial chemicals to 
inhibit the growth of bacteria in the raw Pu-erh tea samples (H. Li et al., 
2018). Cyberlindnera, Uwebraunia, Klebsiella, and Lactobacillus were 
dominant fungal and bacterial genera during the pile-fermentation 
process of primary dark tea (Q. Li et al., 2018), but these genera were 
only barely present in both ripened and raw Pu-erh tea samples, which 
suggested the different microorganisms from different manufacturers 
may serve the same functions for various primary dark tea productions. 
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Fig. 6. Correlations between the contents of chemical compounds and microbiomes. Different bacterial and fungal microorganisms showed a distinct correlation with 
various chemical compounds. Red and orange squares indicated positive correlations between the microbial groups and chemical compounds. (For interpretation of 
the references to color in this figure legend, the reader is referred to the Web version of this article.) 


Bacillus, Rasamsonia, Lichtheimia, and Debaryomyces were mainly re- 
sponsible for flavors production in the late fermentation stage of Pu-erh 
tea manufacturing processes by shotgun metagenomic and metabo- 
lomic analysis (Z. Li et al., 2018). 

Different Microbiomes between the Raw and Ripened Pu-erh 
Teas. The bacterial communities in the ripened Pu-erh tea samples had 
higher diversity than the raw Pu-erh tea samples, while the fungal di- 
versity indices in the ripened Pu-erh tea samples were lower than the 
raw Pu-erh tea samples. The accelerated fermentation process in the 
ripened Pu-erh tea production might have contributed to the microbial 
community differences between the raw and ripened Pu-erh tea sam- 
ples. Results of this study were in general agreement with other studies 


(Ma et al., 2017; Y.; Zhang et al., 2016); while some contradictory re- 
sults (e.g., the discrepancy of alpha-diversity indices, inconsistency of 
bacterial and fungal OTUs numbers (M. Zhao et al., 2015)) require 
standardized methods. Various factors may account for the contra- 
diction and discrepancy among different studies on the microbial 
communities of Pu-erh tea, such as sampling time and size, aging of the 
Pu-erh tea, DNA extraction efficiency, sequencing depth, and coverage, 
production processes, and geography. Furthermore, our findings con- 
firmed that Protobacteria, Firmicutes, and Actinobacteria phyla were 
the core bacteria in the ripened and raw Pu-erh tea (Ma et al., 2017; Y.; 
Zhang et al., 2016; M.; Zhao et al., 2013). 

Pathogens and Beneficial Microorganisms in Pu-erh Teas. 
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Table 3 
Top three microbial predictors for chemical compounds. 


Chemicals Predictor Contribution Portion Rank 
TA Actinobacteria 4289.37 0.25 1 
TA Fungus_Unclassified 2757.51 0.16 2 
TA Cyanobacteria 2334.47 0.14 3 
CAF Actinobacteria 131.14 0.18 1 
CAF Fungus_Unclassified 112.37 0.15 2 
CAF Basidiomycota 109.41 0.15 3 
TPS Actinobacteria 865.92 0.23 1 
TPS Ascomycota 743.33 0.20 2 
TPS Basidiomycota 640.93 0.17 3 
CGA Cyanobacteria 3.67 0.12 1 
CGA Tenericutes 2.91 0.10 2 
CGA Actinobacteria 2.64 0.09 3 
CA Bacteroidetes 164.60 0.15 1 
CA Basidiomycota 133.58 0.13 2 
CA Verrucomicrobia 121.29 0.11 3 
GA Cyanobacteria 78.12 0.12 1 
GA Spirochaetae 66.19 0.10 2 
GA Basidiomycota 59.06 0.09 3 
GCG Gemmatimonadetes 4.94 0.14 1 
GCG Proteobacteria 4.42 0.12 2 
GCG Tenericutes 3.60 0.10 3 
EC Ascomycota 85.27 0.35 1 
EC Basidiomycota 38.07 0.15 2 
EC Actinobacteria 36.78 0.15 3 
Cc Ascomycota 11.61 0.16 1 
Cc Basidiomycota 11.04 0.16 2 
Cc Actinobacteria 8.11 0.11 3 
EGC Bac_Unclassified 43.34 0.16 1 
EGC Ascomycota 23.43 0.08 2 
EGC Proteobacteria 23.04 0.08 3 
ECG Spirochaetae 258.73 0.16 1 
ECG Firmicutes 198.40 0.13 2 
ECG Proteobacteria 160.44 0.10 3 
EGCG Firmicutes 524.24 0.22 1 
EGCG Bacteroidetes 359.53 0.15 2 
EGCG Cyanobacteria 343.02 0.14 3 
GC Firmicutes 9.31 0.10 1 
GC Chlorobi 8.54 0.09 2 
GC Bacteroidetes 7.88 0.09 3 


Another aim is to diagnose the pathogens responsible for secreting 
toxigenic chemicals and identify the probiotics beneficial to human 
health in the Pu-erh teas, which has an impact on the safe production 
and consumption of Pu-erh tea. Aspergillus section Nigri was introduced 
recently, and A. niger fermentation is “generally recognized as safe” 
(GRAS) by the United States Food and Drug Administration. However, 
not all species in the section of Nigri are A. niger sensu stricto. A closer 
screening of Aspergillus and other fungi is required to distinguish the 
beneficial fungi from those mycotoxin-producing, harmful fungi. Haas 
investigated the fungi and mycotoxins from 36 Pu-erh tea samples and 
found that fungi belonging to 19 fungal genera and 31 species could be 
isolated from all Pu-erh tea samples. The most prevalent fungal species 
were Aspergillus acidus, Aspergillus fumigatus, Zygomycetes and 
Penicillium. No aflatoxins and fumonisins were detected in all samples, 
while ochratoxin A was detected in 4 out of the total 36 samples (Haas 
et al., 2013). Pathogens (e.g., Clostridium perfringens) and ochratoxin A 
produced by molds (especially Aspergilli), potentially harmful to the 
human health, were detected in commercial black and green tea sam- 
ples (Carraturo et al., 2018). Despite the detection of harmful micro- 
organisms in the tea samples by the sequencing-based method, species/ 
strain variations may confound the interpretation of biochemical cap- 
abilities of these bacterial and fungal microorganisms in the Pu-erh tea. 
For example, Staphylococcus aureus can produce staphylococcal en- 
terotoxins, causing common food-borne disease, but the Staphylococcus 
genus, the most abundant in the ripened Pu-erh tea samples, were 
different species and may not cause food safety problems (Ma et al., 
2017). Li et al. developed quantitative PCR methods to analyze 
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probiotic bacteria (e.g., Enterococcus, Bacillus, Lactococcus, Bifido- 
bacterium Lactobacillus) during the Pu-erh tea fermentation (Li, Li, Feng, 
& Zhang, 2018). Therefore, improved identification methods and re- 
vised systematics can assist in enhancing microbial safety concerning 
the raw and ripened Pu-erh tea (Houbraken & Samson, 2017). 

In summary, we investigated the microbial and chemical profiles of 
the raw and ripened Pu-erh tea samples and demonstrated the corre- 
lations between microbiomes and chemical constituents in the Pu-erh 
teas. Further studies including cultivation, meta-omics (e.g., meta- 
transcriptomics, metabolomics, proteomics, lipidomics, glycomics) will 
further unravel the complexity of microbial-mediated Pu-erh tea pro- 
duction processes, therefore improving the quality, safety, and benefits 
of Pu-erh teas. 
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